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ABSTRACT 

We provide precise fitting formulae for Faraday conversion and rotation coefficients in 
uniformly magnetized relativistic plasma. The formulae are immediately applicable to 
Rotation Measure and Circular Polarization (CP) production in jets and hot accretion 
flows. We show the recipe and results for arbitrary isotropic particle distributions, in 
particular thermal and power-law. The exact Faraday conversion coefficient is found 
to approach zero with the increasing particle energy. The non-linear corrections of 
Faraday conversion and rotation coefficients are found essential for reliable CP inter- 
pretation of Sgr A*. 

Key words: plasmas - polarization - radiation mechanisms: general - radiative 
transfer - Galaxy: centre 
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1 INTRODUCTION 

The cyclo-synchrotron emission, also called magneto-bremsstrahlung emission, is one of the most important radiative mecha- 
nism in astrophysics. It is believed to produce radio emission in the centers of AGNs and LLAGNs (low luminosity AGNs, such 
as the Galactic Center). A polarized nature of cyclo-synchrotron emission is of increasing interest for radio observers. With the 
help of polarization one can understand the magnetic field structure in radio sources. Basic theory of emissio n and propagation 
of polarized light has been established fe.g jLegg fc Westfoldlll968l : |Pacholczvklll970l : IJones fc Q'Delllll977i ). Particles in cold 
plasma emit cyclotron radiation, which is circularly polarized (CP). When linearly polarized (LP) light propagates through 
cold magnetized plasma, it undergoes Faraday rotation. In turn, relativistic plasma emits synchrotron radiation, which is 
linearly polarized. Light traversing relativistic plasma undergoes both Faraday conversion and Faraday rotation. 

In simple theory the strength of Faraday rotation effect is proportional to \ 2 n e B ■ 51, where A is the photon wavelength. 



is elect ron dens i ty, B is magnetic field vector, and I is the displacement along the line of sight. However, iTrubnikov 



1 195ct ) and lMelrosel (l997al ) 'lave shown that in a general case Faraday rotation of plasma depends also on the Lorentz factor 
of electrons 7. Faraday rotation weakens with the increase of 7 as ln7/7 2 . The electric vector position angle (EVPA) of 
LP light will be preserved better, if the electrons are relativistic. Thus we can infer the intrinsic EVPA of a synchrotron- 
emitting source. The electrons in the vicinity of a black hole Sagittarius A* (Sgr A*) in Galactic Center are often m odeled 
by a relativistic Maxwellian (thermal) distribution with temperature > 10 10 K (e.g. lYuan. Quataert fc Naravanlliooi 'l. When 
Faraday rotation is strong and non-uniform across the beam, then beam depolarization can diminish the resultant LP fraction 



I Bower et al"]l2005l ). The effect of large Lorentz factors on Faraday rotat ion measure near Sgr A* m ust be considered to 



explain the detected LP fluxes in sub-millimet er to near- infrared bands ( Quataert fc Gruzinov 20ld ). Previous work (e.g. 
Ginzburg I Il964l ; ISazonol Il969l : [Melrose! Il997bl ) also highlighted Faraday conversion, or generalized Far aday rotation. This 
quantity describes the interconversion of linearly and circularly polarized light. It is normally expected ( Homan et al. ,2009) 
to convert emitted LP light into CP light during propagatio n in m agnetized medium. This is the likely cause of o bserved high 
CP fraction of Sgr A* spectrum in radio band I Bower et al. 20021 ) and sub-millimeter band ( Munoz et al. 2011 ). 
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The strength of Faraday conversion was found proportional to A 3 n e £> 2 sin 2 951, where 9 is the angle between B and 81 
(or k, see in Figure [l}. Additional suggested proportionality to electron temperature T e makes Faraday conversion reach very 
large values in r elativistic plasmas. However, this proportionality ceases at very high T e and Faraday conversion measure 
approaches zero Jshcherba cov[ 2008h . A detectable CP fraction can be generated near Sgr A* in sub-mi llimeter band (e.g 
Balla ntvne. Qzel fc Psaltisl l2007l )7 but precise treatment of Faraday rotation and conversion is essential (jShcherbakov et al. 



2010h . A non-zero CP fraction is already detected with SMA (Munoz et al., 2011, ApJ, submitted) at 230 GHz and 345 GHz 



This by itself points in the direction of a v ery hot radiatively i neffici ent accretion flow (RIAF). The origin of CP near Sgr A* 
and jets was quantified bv lBeckert |2003h ; IShcherbakov et all (|2010h . Yet more detailed and accurate calculations are needed 
to q uantify the circular polarizations for the Galactic Center (GC) supermassive black hole and other radio source. 

Huang et al. (2009) was the first to incorporate Faraday rotati on and conversion withi n gene ral relativistic (GR) polarized 



radiative transfer framework, though with some approximations. IShcherbakov fc Huang] ( 201]J , hereafter Paper I) provided 
a method to accurately calculate cyclo-synchrotron absorption and Faraday conversion/rotation for electrons in isotropic 
thermal distribution and outlined the exact procedure for polarized radiative transfer i n GR. Precise Faraday rotation and 
conversion coefficients were computed earlier for thermal plasmas in Shcherbakov 1 20081 ). An important cornerstone in com- 
puting propagation effects is linear approximation, in which only the first non-zero terms in series expansion in Q.q/lo ratio 
ar e taken for corres pondent quantities (e.g. formula l47)l . It was also derived therein and found c onsistent with result provided 
bv lMelrosel (|l997bh . The precise values of Faraday conversion coefficient in IShcherbakov! {2008) and Paper I match the linear 
approximation from cold to weakly-relativistic regimes of thermal plasma. In relativistic regime Faraday conversion largely 
deviated from the linear approximation, because it breaks at finite ratio of £Iq/uj. Here fio = eB/(m e c) is the cyclotron fre- 
quency and uj is the radiation frequency. Due to the lack of full investigation of various electron distributions, the important 
question was left unanswered: is thermal distribution special or such behavior of Faraday conversion is generic? 

In present paper we expand the computations of Faraday conversion and rotation coefficients to non-thermal particle 
distributions. We compute the absorption coefficients as well with the same unique method. We find solutions of the wave 
equation and natural modes from cold limit to ultra-relativistic limit. Our formulae are precise at all reasonable particle 7's. 
We find a large discrepancy, if the linear approximations to Faraday conversion and rotation are used, thus justifying the need 
for precise computations. To be practical we adopt 5-function energy distribution of electrons and provide the fitting formulae, 
which can then be integrated over any isotropic distribution of particles. We also provide public code in Mathematica 8 to 
numerically compute the integrals. The paper is organized as follows. We derive dielectric tensor and dispersion relations for 
uniformly magnetized relativistic plasma with isotropic monoenergetic particle distribution in §2. The properties of natural 
modes are investigated in §3. In §4 we provide simplified formulae and generalize to arbitrary electron energy distributions, on 
its polarization prediction is described In §5 we show that the polarized spectrum changes a lot, when linear approximations 
are used for Sgr A* modeling. 



2 RESPONSE TENSOR AND EIGENMODES OF UNIFORMLY MAGNETIZED RELATIVISTIC 
PLASMAS WITH 5-FUNCTION ENERGY DISTRIBUTION 

2.1 Geometry and definitions 

Let us define a coordinate system in three-dimensional flat space with ei along the major axis of the synchrotron radiation 
ellipse, e.2 along the minor axis, and towards the observer (see Figure[l|. Vector ei is perpendicular to B. This coordinate 
system is right-handed, i.e., the observer finds a counter-clockwise rotation from ei to e 2 . We define Minkowski space-time 
with basis eJJ, so that = (—1,0,0,0) and = (0, Bj),j = 1,2,3. Unlike in Paper I, here we perform all derivations in 
spatial basis (ei, e2, 63), instead of (ei, §2, 63). In latter basis ei = ei and 63 || B. We set the speed of light to unity c = 1. 
Normalized vectors of electric field in a transverse wave are 

Ei = ei e l{kx) and E 2 = e 2 e l{kx) , (1) 

where fc M = lo{— 1, 0, 0, 1) T is the covariant photon momentum, are four-coordinates, and kX represents the inner product 
kX — fc M A M . The projections of electric field E of an arbitrary wave along these unit vectors have complex amplitudes A 1 ' 2 , 
or real amplitudes A 1 ' 2 and phases 6 ' 2 , so that 

E = E 1 + E 2 = A 1 E 1 +A 2 E 2 = A 1 ^^ Ei + A 2 e u '' '' E 2 . (2) 
The tensor of intensity is defined as 

Tij =<E,- conj[E 3 } >, (3) 



where conj[...] stands for complex conjugate and <> represents the average over the wave ensemble. 
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Figure 1. Geometry of the problem. Vector B represents uniform magnetic field. The transverse plane wave travels along k and has 
electric field E in (ei,e2) plane. 



Radiative transfer in uniform medium is described by the equation I Sazonov 19691 ) 



lo d 
2tt ' ds lJ 
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(4) 



where s is the distance along the ray, ey is the tensor of spontaneous emission and a 13 is the tensor of wave propagation, or 



the response tensor. The equation ([4]) can be rewritten in a more familiar form IjShcherbakov fc Huanal201ll ) 
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with the polarization vector 
S= (I,Q,U,Vf 



(5) 



(6) 



being the vector of Stokes parameters. The intensities S and integrated polarized fluxes can be directly observed. Here 
Ei,EQ,ev are the emission coefficients, 

rji = Im(a 22 + a 11 )/^, 
tiq = Im(a 1 — a 22 )/^, 

Vv = 2Re(a 12 )/u, (7) 
are the absorption coefficients, 
pv = 2Im(a 12 )/^, 

PQ = Re(a 22 - a 11 )/^, (8) 



where v = u/(2ii) is the frequency. In the following we will concentrate mainly on Faraday rotation coefficient pv and Faraday 
conversion coefficient pQ, which are generally called propagation coefficients. They directly influence the observed polarized 
fluxes. 



2.2 Response tensor and dispersion relations 

We start with the formula for the response tensor for isotropic electrons 

cT{k) = -— I d 3 p^ZW + ^ r d££(-£) 
m e c J d7 m e c J Q 

where W = (1, 0, 0, 0), W is 4-velocity of electrons in observer's Minkowskian frame, p is dimensionless 3-momentum defined 
as p = y/y 2 — 1, and f(-y) is the energy distribution function of electrons. It is normalized to the number density of electrons 

n e as 



d 2 



d 3 / _Id/W e -««)u+sw 
7 d7 



(9) 
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f(j)d p = n e 



(10) 



Tensor f M "(£) describes how the velocity of electrons changes with proper time 1Z(£)U represents the differe nce of the phase 
(kX) of the electron, and S u i s an au xiliary var iable. The formula © coincides with Eq. (2.3.11), Chapter II. iMelrosel l|201(]h 
and with Eq. (19) in iMelrosd l|1997al ). except for £ M "(£), because we defined a different coordinate system. Eq. (35) in Paper 
I offers a similar expression derived in the 3-dimension space with basis (ei,e2,e3). The derivation of Eq. [9] and the related 
definitions can be found in Appendices A & B. 

The 4-vectors in expression ([9]) can be split into temporal and spatial parts as 



72. M = (-u£,Rf), W'* = (7,u i ), S„ = (s ,s j ). 



(11) 



Then the phase becomes —ilZ(£)U + SU — ■ 7 — tKjU 3 + S07 + SjU 3 = — JSo)7 — «(Rj + JSj)u J . The momentum integral 
in the response tensor is 

d » p f 1 d/(7) el( ^_ lso)7 _, (R . +1B . )u j 
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The absolute values are taken as R + is\ = -J (R fc + is k )(Rk + iSk) and |u| = \/u fc Ufc. Note that for any distribution with 
/(7min) ■ j4(7min) — > and /(7 max ) ■ A(7 max ) — » 0, the second term in last expression of Eq. (fi"2"j) vanishes. 

As the first step of calculation for arbitrary distribution of electrons, we use a 5-function as the distribution 



/( 7 ) =5(7-70). 

Then 7(^,5) becomes /(70; £,,S) as 
/(7o;£,S) = e '^-"o)~/o 



sm(R + zsp ) ,7o nn , \~ \ 
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where 

Po = Vol - 1 



(13) 
(14) 
(15) 



is the dimensionless momentum. 

Now we apply the differential operator d 2 / (dS^dSa) to 7 (70; £, 5), set 5 M = 0, and only choose (p, v) — (1,2) to isolate 
the transverse wave component. The final expression of 2 x 2 response tensor is 
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t«(0 
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R ! (£)R^) 
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The differentiation of -A (7; £, 5) yields a response tensor boundary term as 
47re 2 



- cos 9 (sin(fi £) - fi £) (1 - cos(OoO) 
cos 2 (sin(O £) - ^oO 2 
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The expression ()17|) is needed for distributions confined by cut-off Lorentz factors. 
The dielectric tensor 



5 13 + 



Aire 



leads to the wave equation 



k 2 c 2 



(17) 



(18) 



(19) 



(20) 



The corresponding wave dispersion relation is 
k\<? = LO 2 + 2-kc a 11 + a 22 ± \J (a 11 — a 22 ) 2 + 4a 12 a 2 

where a 21 = — a 12 due to Onsager principle ( Landau fc Lifshitsl 1980h (p. 273). Note that k\ are real, when ol 3 only has 
the Hermitian part, while k\ become complex, if of 3 has both the Hermitian and the anti-Hermitian parts. The approximate 
relation k\.(? w uJ 2 + (27rc)Re(a 11 + a 22 ) helps to det ermine if the waves can propagate through medium. For example, the 
number density n e < 10 7 was estimated near Sgr A* ( Yuan. Quataert fc Naravan 20031 ) . so that k± are dominated by their 
real parts for uj > 100MHz and the waves can propagate. 



NUMERICAL CALCULATION OF RESPONSE TENSOR AND THE ELLIPTICAL NATURE OF 
NATURAL MODES 



3.1 Numerical calculation of response tensor 

We substitute the components of the response tensor from Eq. (|16|l and perform integration over £ in complex plane in 
Mathematica 8. The source code in Mathematica 8 can be found at http://astroman.org/Faraday_conversion/ We then 
find propagation coefficients pv and pq and absorption coefficients rji , t)q, and rjv according to relations (|7I8[) . We analogously 
compute the boundary terms by substituting the components of «g from Eq. (|17|l . Just like in Paper I, we do not perform the 
integration over £ along the real axis. To accelerate the convergence we integrate in a complex plane along the ray originating 
at £ = at a positive angle ij) G (0,7r/2) to the real axis. Angle ip needs to be small enough in order to avoid crossing the 
branch points of 01 3 . These branch points are produced by zeros of R(£)- We integrate the full complex expressions to find 
the response tensor a 13 and the boundary term. If 70 is small, then computations of anti-Hermitian parts involve substantial 
cancelations with the values of integrals close to zero. Thus it is hard to reach good accuracy for absorptivity calculations 
with the method chosen, whereas the correspondent Hermitian parts easily converge. When 70 is larger than ~ 10, the values 
of Hermitian-related part and anti-Hermitian-related part become comparable and all integrals converge. 

We choose the fiducial model with Q,q/lu = 10~ 4 and 9 — 45° and plot on Figure __| propagation coefficients and absorption 
coefficients as functions of the Lorentz factor. On the left panel, 2Im(a 12 ) oc pv and — Re(a n — a 22 ) oc pq, multiplied by 
(m e c)/ (47re 2 ), are shown in red and purple dotted lines, respectively. They both monotonically increase as 70 increases from 1 
to ~ 60. As 70 increases further, the profile of pv becomes flatter, while pq reaches its peak and decreases to negative values. 

On the right panel, lm(a 11 — a 22 ) oc tjq and 2Re(a 12 ) oc rjv, divided by (47re 2 /m e c), are shown in green and blue dotted 
lines, respectively. We also show ( u;/27r)r?o and (uj/2ir)r]' v as long-dashed and dashed black lines, respectively. Here r/'g and 



rjv corresponds to the integrals in lSazonov 
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(21) 



1 The sign of rjv is opposite to that in lSazo nov ( 1969|), if the current IAU/IEEE definition of the sign of circular polarization is followed. 
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Figure 2. Absorption and propagation coefficients for the fiducial model with parameters Qq/uj = 10~ 4 and 6 = 45° for the monoenergetic 
electron distribution. Left: 2Im(o 12 ) oc py (Faraday rotation, in red) and — Re(o 11 — a 22 ) oc pg (Faraday conversion, in purple). Right: 
lm(a 11 — a 22 ) oc t)q (absorption of linearly-po larized w a ves, g reen) and 2Re(o 12 ) oc r]y (absorption of circularly-polarized waves, blue). 
Approximate absorption coefficients based on ISazonovl l|l969t ). in particular, (aj/27r)r)g and (u)/2n)r]' v , are shown in long-dashed and 
dashed black lines, respectively. 



where vb is the cyclotron frequency, f 7 = (3eB sin 8j 2 ) / (4irm e c) is the characteristic frequency, and K a is the modified Bessel 
function of the second kind of order a. We integrate by parts to deal with the differential of the 5-function. Note that our 
absorptivities deviate a lot from the correct values, when 70 is small. As mentioned earlier, this is due to large cancelations 
of the parts of the integral, so that absorptivities cannot easily converge. When 70 > 30, they coincide with the approximate 
expressions. 

Despite inaccuracy at low Lorentz factors, the calculation clearly shows an important property of plasma absorption. The 
absorption in Q-component is smaller, compared to that in V-component r)Q < \r]v\ at low Lorentz factors, while t)q > \r)v\ 
at high Lorentz factors. This shows that the radiation mechanism changes from CP-dominated cyclotron to LP-dominated 
synchrotron as the Lorentz factor increases. The traditional approximations do not exhibit this property, because they assume 
high Lorentz factors. 



3.2 The axial radios and natural modes 

The corresponding eigenvectors from the wave equation are (T + ,t) T and (T~ ,i) T , respectively, where 



^ + y /(a 11 -a 22 ) 2 + 4a 12 a 2 



2ia 12 

~_ a 11 - a 22 - ^/(a 11 - a 22 ) 2 + 4a 12 a 21 



(22) 



2jq 12 

These T + and T~ are complex axial ratios of the transverse wave. They obey the relation 
f + f~ = -1. (23) 
The polar decomposition into real amplitudes T 0) and phases ip reads 

f+ = T + e ^ + and f~ = T~e ltp ~, (24) 
so that 

T + T~ = 1 and tp + + <p~ = -tt. (25) 
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These two wave eigenvalues define two natural wave modes as 



T+E! + iE 2 TV V+ Ei + iE 2 



E+ = 



1 + T+ ■ conj[T^ 



VI + (T+) 2 



~ T~Ei+iE 2 T~e llp E!+iE 2 , oc , 
E = = , ,„ (26) 



1 + T~ ■ conj[T-] \/l + (T f 



In practise, tp + < —ir/2 and tp > —ir/2. Thus, the electric vector E+ rotates counter-clockwise and the electric vector E- 
clockwise, as seen by the observer. Note that 

#, . r " . e <f+-f-) + 1 ( =0, iff <p+=-n 

E + ■ conj\E-} = : — < , ' , (27) 

+ Jl 1 T++T- \ t^O, otherwise v ; 

That is E+ and E- are not perpendicular to each other, unless the anti-Hermitian part, absorption, can be neglected. 



3.3 Properties of radiation in natural modes 

Let us decompose the natural modes along Ei and E 2 as 

E ± = E] and £ ± = Bo (28) 

a/i + (t±) 2 v/iT(r±f 

We define the Stokes parameters for natural modes as 



7 ± = <Ef ■conj[Ef]> + <Et-conj[E 2 t ]>=(A ± y 



= <Ef- conj[Et] >-<Et conj[Ef] >= (A*)' 



i + (T±y 



Uf = <Ef- conj[Ef] > + <Ef- conj[ E f] >= (A*) 2 ■ , 

V* = x(<Et- conj[Ef] >-<Ef. canj[Ef] >) = (A^ 2 ■ 2 ^ ^ . (29) 

Note that .Xo 7^ -X^~ + Xq,(Xq = Jo, <3o, Uo, Vo), since < E+ ■ conj[E-\ >^ unless <p + — —iv. These Stokes parameters 
correspond to elliptically polarized radiation with ellipticity j3 

fl± 1 • -1^ 1 ■ -i/2T±cos^±\ 

a = — sm — t = - sin rr (30 

2 j± 2 ^ 1 + (r±) 2 y V ' 

and the electric vector position angle (EVPA) 

X* = itan-4 = itan-^ 2r± f (31) 
X 2 Q± 2 ^ (T±) 2 -1 J K ' 

The relations between the axial ratios and the Stokes parameters for each mode are 



jRe(T + )Re(T")|. (32) 



Kp) V (Oo + )+(^)7 fe) V (Qo)+WV VIT++T- 

We define a special quantity 

Z ± = ^-^, (33) 

which effectively measures how well the normal modes can be described, if we ignore absorption. We show amplitudes T ± (^ 0), 
phases tp , ellipticities /3 , and EVPA x on Figure [3] All quantities with superscript ( + ) are drawn in solid lines and those 
with (~) are dashed. 
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Figure 3. Properties of two natural modes for the monoenergetic particle distribution. Amplitudes (black) and phases ip^ (grey) of 
complex axial ratio, ellipticity 0^ (~ ±7r/4 corresponding to circular polarized waves and ~ to linear ones), EVPA x > and = Z~ 
(eg. 1331) arc shown in four panels, respectively. Lines for (+)-mode are solid and for (~)-mode are dashed. 



3. 3. 1 Cold plasma limit 

In cold plasma the anti-Hermitian part of the response tensor can be neglected compared to the Hermitian part. In this case, 
(p + — — 7r, (p ^ — (see grey lines on the top left panel of Figure[3|, Uq = Uq = 0, and \py\ >> \pq\. The axial ratios are 



T 



Pv/PQ 



-T 



1 + V1 + (Pv/pq) 2 

Pv/PQ 



(34) 



They are consistent with the Eq.(4.6) in lMelrosd |l997bj). The corresponding ellipticity /3 and EVPAs are 
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+ " 
x = -3, x 



= 0. 



(35) 



Two natural modes are circularly polarized. They are orthogonal with major and minor axes aligned with B\ and respec- 
tively. The relation between the Stokes parameters and the axial ratios becomes 
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,± 
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-Rcp 
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Qo 



T+ - T~ 



pq' 



1. 



(36) 



This means the ratio of circular to linear radiation intensities in cold plasma eigenmodes equals the ratio of Faraday rotation 
to Faraday conversion coefficients. The total emissio n is dispersionless i n a sense that the te rm (p x p) in transfer of polarized 
radiation in lLandi Deel'Innocenti fc Landolfil (J200J) vanishes, also see iHuang et all (|2009h for similar results with different 
definitions of the axial ratios. As Faraday rotation is much stronger than Faraday conversion, then the circular polarized inten- 
sity is much larger than the linear polarized intensity in eigenmodes. These are the well-known properties of 'cyclotron'/cold 
plasma regime. 

In the limit lu 2 >> a?" , the refractive indices of two modes are 



k±c 



2nc 



1 2nc, 



n± = ^ = Wl + ^(« 11 +^±^pv)^l + ^(a 11 +a 22 ±^pv), 
lu y ur 2tv 2 lj^ 2iy 

which corresponds to standard rotation by Of of EVPA plane, Faraday rotation effect, as 

dO 



QPf LU . 

— - = — [n+ — ri- 
dl 2c y + 



2 PV 



(37) 



(38) 



3.3.2 Ultrarelativistic plasma limit 

In plasma with very high Lorentz factors pQ changes its sign to the negative. Two natural modes also change their polarization 
property from circular to linear. We find an interesting result that the quantity Z is approaching another constant. This 
means the proportion of linear to circular radiation in each eigenmode can still be easily measured by the proportion of 
Faraday conversion to rotation. 

As it appears t\q « —2pQ in this limit, and \r^v\ « \pv\ « \pq\- So the complex axial ratios become 



T= 



PQ 



PQ 



1 PQ 



PQ J \PQ J 



PV 



[(£)+»] 



pv 




We derive 




T + -> 


PQ ( 




Pv 




_PV 


PQ 


± -> 


0, X 


and 





l + 2t± y/(l + 2if + (pv/pq) 2 



PQ_ 

PV 



T + ->• 00, lp + ^ -it + tan x (2)^-2, 
.2i), T^^O, if' = -7T - ip + ->• -tan _1 (2) 

_ 7T 

x -> -3, 



CP 



CP _ 
'LP 

-0.2 



\T+ +T- 



Re(T+)Re(T")| 



4-47- I pq/pv IV Pv 



2^-0.1^ 



PQ 



(39) 



(40) 



-0.2-^ 

PQ 



(41) 



Note that the total intensity cannot be calculated by simply adding intensities in two modes, because they are not orthogonal. 



3.3.3 Intermediate regime 

In plasma with intermediate Lorentz factors (about 70 ~ 10 1-2 for the fiducial model with flo/w = 10 -4 ), the properties of 
the natural modes change gradually from those in cold limit to those in ultra-relativistic limit. 

As 70 increases, T + decreases to 0, while approaches 00, ellipticities \f3 \ decrease to 0, EVPAs x ± deviate from -90° 
(or 0°). At a special value of 70 the Faraday conversion coefficient — Re(a 11 — a 22 ) changes its sign, as shown in Figure [2] 
The Lorentz factor for the fiducial model is 70 ~ 100. Flips of T , f) , and x between two modes help to preserve the 
handedness of the modes. That is the wave with E+ always has a counter-clockwise rotation, while the wave described by 
E- has a clockwise one. The value of Z decreases in the intermediate regime from 1 to ~ —0.2. 
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4 PROPAGATION AND ABSORPTION COEFFICIENTS FOR ISOTROPIC PLASMAS WITH 
ARBITRARY ENERGY DISTRIBUTIONS 

4.1 Integration over particle distribution function 

The final form of the response tensor is related to I(£,S) as shown by Eq. (1121) . With the aid of 
= y"/(7)/(7o;C,5)d 7 , 

one has 



(42) 



2 /*°° 

m e c J 



d 2 ( d 2 



dSudS, 



,f( 7 ) ■ a IJ 0, 7 )d 7 - /( 7m ax) • QB J ( fc '7max) +/( 7 min) • QB( fc '7min). 



(43) 



The isotropic distribution function /( 7 ) is normalized as J /(7)d 3 p = n e . In this section we compute the propagation and 
absorption coefficients for thermal and power-law particle distributions. 



4-1.1 Thermal distribution 



The thermal distribution is 



/(7) = 



-7/e e 



1 < 7 < +oo, 



where O e = kBT e /(m e c 2 ) is the dimensionless particle temperature. It is normalized to number density of electrons n e as 



(44) 



/( 7 )d 3 p = 4-7T J jy/ 7 2 — l/( 7 )d 7 = n e . 

In this case, /( 7max — > +oo) — >• 0, so that the response tensor becomes 
r />+oo 

qI3 W = 77^7?wT / e-^ e --a«(* )7 )d7 + e- 1 / e -.ag(*,l) 

47T(s) e A2(fc) e ) [./I 



47re e A- 2 (e- 1 ) 



(45) 



- 7/e =-a y '(fc, 7 )d T , (46) 



since Qg(fc, 1) — ¥ 0. The propagation and absorption coefficients for n e — 1 are shown in thick dotted lines on Figure [J] 



On the left panel we sho w in dashed grey lines the linear approximations to propagation coefficients elaborated in Melrose! 
1 1997bl l; IShcherbakovl (120081 ) for n e = 1. The related formulae ar«Q 



PQ, Hn = 



PV.lin 



Ki(e; 



27re 2 no 



4ire 2 n KoiQ- 1 ) 
m e c-iu 2 K 2 (Qe 1 ) 



■ sin 2 9 + 60 e sin 2 I 



cos 0. 



(47) 



We scale 2Im(a 12 ) and (o;/27r)pv,ii n by a factor of 300 for a better layout together with Re(a 11 — a 22 ) and (u/2w)pQ l u n . Note 
that py,iin is a good approximation for Faraday rotation coefficient at any temperature, while PQ,n n is a good approximation 
for Faraday conversion coefficient only at low temperatures, not at high temperatures. The peak of Faraday conversion for 
monoenergetic particles leads to a similar peak at about 10 11 K for thermal distribution. Faraday conversion is much lower than 
the linear approximation predicts, if the temperature rises. We will show in the last section that the linear approximations 
lead to wrong predictions of circular polarization from Sgr A*, the difference being a factor of several. 

Black long-dashed and dashed lines on the right panel show (cj/2TT)riQ and (u)/2ir)r)' v , for n e = 1, respectively. The 
absorption coefficients r/'g and rj' v are calculated according to Eq. (|21[) by substituting the thermal distribution function. 
Similar to the case of monoenergetic particles, our thermal calculations match the approximations well at high temperatures, 
while they become inaccurate in low temperatures. In practice one can just adopt the simple traditional approximations for 
synchrotron absorption coefficients. 

We show T and phases, f3 , x , an d Z on Figure [5] In general, they are similar to those for <5-function distribution 
discussed in §[331 Note that EVPA (x) deviates by as much as ~ 20° from 0° (or -90°) at T e w 10 11 K. 



2 Note, that Faraday conversion coefficient in the "linear" regime is actually proportional to [Q,q/uj) 2 . 
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Figure 4. Propagat ion and abs orption coefficients for thermal energy distribution. Left: Faraday rotation and conversion. Linear ap- 
proximations bv lMelrosel jl997bh are shown in black dashed lines. Right: Absorption coefficients. Traditional approximations are shown 
in black solid lines. 



4- 1.2 Power-law distribution 

Number density per unit Lorentz factor is 
b 



N(j) 



i mm 



_l - b 

/max 



7min < 7 < 7n 



,b> 1 



(48) 



for power-law particle distribution, where 6 is the energy-spectral index, set greater than 1 as examples. Thus the distribution 
function is 

fil) = - A . i = 1 1=6 1=6 • J 1 ==, 7min<7<7max. (49) 



We set 7max = +oo in computations below for b > 1. Then /( 7n 
n e (6-l) 



+00) — > 0, and the response tensor becomes 



a 3 {k) 



a l:l (fc,7)d7 + 



-b 
7min 



V ' mm 



«B ( fc > 7min) 



(50) 



We show the results of num erical integration for b = 2.5 on Figure [6] On the left panel, we also show approximations for 
propagation coefficients given bv lSazonovl (jl96^ 1 for n e = 1 in dashed grey lines. The related formulae are 



PQ.appr = 8.5 x 10 



1.7 x 10"" 



(6-2) 
In 7 m 



fio sin 6j, 
(6-1 



2 

min 



(b-2)/2 



n 



(6 + 1)7 



5+1 
m i n 



7 n 



- 1 

2tt 



(6-1) (Via 



(b+2)/2 



2tt 

1 



(51) 



(Similar formulae can be also found in I Jones fc O'Delll (|l977T ).) 

We scale 2Im(« 12 ) and (w/27r)/Jv ia ppr by a factor of 150 for a better layout. Here pv>p P r is a good approximation of 
Faraday rotation coefficient only for large 7mm- It significantly underestimates Faraday rotation, if the cut-off Lorentz factor 
is low. The approximation pQ, apP r works well for Faraday conversion coefficient, if 7 m in < 100. It accurately describes the 
exact behavior including the peak. 

On the right panel of Figure [6] long-dashed and dashed black lines show {ijj/2-k)t)'q and (to /2ir)r)' v respectively for n e = 1. 
The approximate absorption coefficients rjq and r) v are computed based on Eq.[2T]by substituting the power-law distribution 
function. In this case our calculation of In^a 11 — a 22 ) and 2Re(a 12 ) match their traditional approximations well for all 7 m i n . 
This is because the particles with high 7 play a large role in power-law distribution compared to the thermal, so that the 
inaccuracy at low 7's is concealed. Similar to the case of thermal distribution, one can adopt simple traditional approximations 
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Figure 5. Properties of two natural modes in plasma with thermal energy distribution.. Amplitudes (black) and phases tp^ 1 (grey) 
of complex axial ratio, ellipticity /3 , EVPA x i an d Z + = Z~ (eg. 133ft are shown on four panels, respectively. 



for synchrotron absorption. We also show approximations for absorption coefficients given by ISazonov (1969) for n e = 1, in 
long-dashed and dashed grey lines. The corresponding expressions are 

, lvln -4 , n ^( 3b + 2 \^( 3b+W \ (6-1) / a, . y b+2)/2 ^ 
%, appI = 3.1x10 .(b + ^r^^-jr^^-j-j^-^-smej ■ -, 

' mm 

|±lr («£!) r (»±!1) . <» + ■ £ . 



^/V.appr 



' mm 



7 n 

These approximations are good for 7 m in < 100, while they overestimate the absorption at larger 7 m i n . 

We show T and phases, /3 , x , and Z in black and grey on Figure [7] for power-law particle distribution with b = 2.5. 
The same quantities for power-law distribution with b — 1.5 are shown in green and cyan for comparison. The plots are similar 
to those for monoenergetic distribution or thermal distribution. The curves for 6 = 1.5 and b = 2.5 almost coincide for high 
7min- At low 7 m i n , however, x deviate more from —90° (or 0°) at b = 1.5 compared to the case with b = 2.5, and Z ± 7^ 1 for 
b — 1.5. This is because the particles with high Lorentz factors affect the results more for a lower spectral index. 
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Figure 6. Propagation and absor ption co efficients for power-law energy distribution with index b = 2.5. Left: Faraday rotation and con- 
version. Linear approximations by [Melrose! <|l997bh are shown in grey dashed li nes. Right : Abso rption coefficients. Traditional synchrotron 
approximations are shown in black solid lines, whereas the approximations bv lSazonovl are shown in grey dashed lines. 



4.2 New approximate formulae for Faraday conversion/rotation coefficients 

We have shown that the traditional synchrotron approximations to the absorption coefficients are accurate and practical. On 
the contrary, si mple linear approximations of the propagation coefficients have large errors at high elect ron energies. Although 
Sazonov ( 1969) also provided the integral expressions for the propagation coefficients in their Eq. (2.3). |jones fc O'Delll (|l977h 



in their Eq. (C16) as well, the Faraday rotation (/' r ' therein) and especially Faraday conversion (/i (r) therein) are inaccurate 
for high energies of electrons. 

We devise a new set of approximate formulae for the complex response tensor in a plasma with monoenergetic particle 
distribution. The goal is to provide simple relations for accurate evaluation of propagation coefficients. We provide the 
numerical code in Mathematica 8 for the full evaluation at http://astroman.org/Faraday_conversion/, but we encourage 
the readers to use the simplified formulae for practical applications. These simplified formulae for Faraday rotation and 
Faraday conversion are computed for plasma with 5-function energy distribution. They can be easily integr ated over the 



arbitr ary energy distribution. Similar approximate formulae for thermal particle distributions were computed in IShcherbakov 



(2008 



4-2.1 Computations for monoenergetic particle distribution 



It is non-trivial that good approximations to Faraday rotation/conversion coefficients exist in a three-dimensional parameter 
space of flo/ui, 9, and 7. However, we manage to find the formulae accurate to within 10% at most reasonable combinations 
of parameters. 

We define an auxiliary quantity 



X A 



^2 sin 6» fi 



10- 



and introduce four new expressions: Hx, Hb, <?x, and gs to approximate — Re(a n 
2Im(ag 2 ), respectively. The formulae are 



Hx(-yo) 



9.29 • 10- 9 ■ y/1 - 7^ • (Aa7o) 3 



-0.000203(Xa7o) 04343 
+0.002 exp 



(In X A 7Q -4.2137) 
0.5429 



0.0013 cos [0.5646 ln(X A 7o 
- 0.00083 exp 



4.03] 



(In X A 7Q -4.2137] 
0.2121 



(53) 

a 22 ), -Re(a^ - a| 2 ), 2Im(a 12 ), and 
Xa7o < 40, 

Xa-jo > 40, 
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Figure 7. Properties of two natural modes for the power-law particle distribution. Amplitudes (black) and phases ip^ (grey) of 
complex axial ratio, ellipticity /8 , EVPA x > an d ^ + = Z~ ( ec l- I33H are shown in four panels, respectively, for b = 2.5. Lines for 
(+)-mode are solid and for ( — )-mode are dashed. The same quantities for the power-law distribution with b = 1.5 arc shown in green 
and cyan. 



#b(7o) = < 



4.67-10- 9 -(l- 7 (7 1 )* -(Aa 7 o) 3 - 84 , 



0.864 - 0.2082 (lnX A 7o) + 0.0175 (lnX A 7o) - 0.000626 (lnXv7o) 



X A jo < 40, 



+1.0175- 10~ 5 (lnX A 7o) 8 - 7 
-0.01 exp 



(In A' A 70-4. 0755) 
0.0763 



px(7o) 

SB (70) 



1 - 0.4 exp 



(In X A7o - 9.21) 



11.93 



0.05 exp 



Hr 8 (lnX A7 o) 10 
(lnX A70 - 5.76) 2 



1.33 



+ 0.075 exp 



X A7o > 40, 

(lnAVyo -4.03) 2 



0.65 



1 - 0.0045(X A 7o) 



(54) 
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Figure 8. Accuracy checking for simplified formulae for Faraday conversion pq and Faraday rotation py. See details in text. 



The Faraday conversion/rotation coefficients for monoenergetic isotropic particle distribution with Lorentz factor 70 are 
calculated as 



PQ (7o) 
PQ,b(7o) 

Pv(7o) 
Pv,b(7o) 



2tt 

0J 

2tt 

UJ 



Re(a 11 
Re(Qe 1 



22n 



8ttV 
m e coj 
8ttV 



X A -#x(7o), 
#b(7o), 



2 . ^l m (a 12 ) 
2--Im(^ 2 ) 



8Tv 2 e 2 n cos 9 
m e cuj 2 

87rVfi cosd 



In 



1 +P0/70 



70 • In 



1 — po/70 

1 +P0/70 



5x(7o), 



1 -po/70 

Those for arbitrary particle distribution 7(7) are calculated as 



2p 



• SB (70) 



(55) 



PQ = 



Pv 



fil) ■ PQ (7)d7-/(7n 
/(7) -Pv(7)d7-/(7n 



PQ,B(7max) + ,f(7min) ' /OQ.B^min), 
' PV,B(7max) + /(7min) ■ pV,B(7min)- 



(56) 



The above formulae are accurate, with errors of several percent for ^-distribution in general. As a test, we choose f("f) 
in thermal distribution and calculate (cj/2iy)pq and (cj/27r)p v by them for a wide range of temperature. We then calculate 
(u}/2tt)pq and (uu/2ir)py by contour integrate formulae for thermal distribution, which is discussed in the Sec. 4. 2. 2. The errors 
of pq are defined as \(uj/2tt)pq — (uj/2tv)pq\/\(uj/2ty)pq\, and those of pv, similarly, as j(tj/27r)pv — (w/27r)/9 v |/|(w/27r)p v |. 
I.e., error of 1 means the coefficient is accurate with a factor of 2. We choose nine couples of parameters (-^f-,6 1 ), such as 



(10" 2 ,l o ), (Hn 2 ^ ), (10- 2 ,89°), (10" 4 ,1°), (10- 4 ,45°), (10- 4 ,89°), (HT S ,1°), (Hn 8 ,45°), and (10" 



as examples 



represented by different symbols (cross, square, and triangle) and different color (blue, green, and red), respectively. As shown 
on Figure [8] the errors are less than 1 (100%) in general. For high temperatures at which synchrotron emission is effective 
(> 10 10 K), the errors are as good as within 0.1 (10%). 



4-2.2 Results for thermal distribution 



We have shown in § 14.1.11 that one can calculate the response tensor for monoenergetic particle first by integrating over the 
proper time (or phase), then integrate the result over the arbitrary particle distribution. For thermal distribution, one can 
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change the integration order and analytically integrate over 7. Thus, the response tensor can be simplified to 



a J (k) = 
where 
K'(Q = 



in e e 



m e cQ- 2 K 2 (Ql 1 ) 



(57) 



9~ 2 - 2ie~ 1 uj^ + 



()-' 

"0 



-[2 - 2cos(n £) - n 2 ^ 2 ] 



and t l J (Q and T lJ (£) a re gi ven in § 12.21 Sec the corresponding derivation in Appendix C, which closely follows iTrubnikov 
1 1958U Melrose] (|l997al l. and lSwansonl ~(f2003). The integrals for thermal a %i (k) in Eq. (f5T|) are easier to compute numerically. 
However, such simplification can on ly be done for thermal distribution. The detailed computations and discussion of Eg. 1571 
can be found in Shcherbakov (2008). Good fittings for thermal distribution, accurate within 10% except with large Qo/u and 
large 6, are also provided therein. We will provide their expressions for electrons in a present notatiorjj: 



PV,th = 



7; ^yisin^lO 3 ^) 



PQ.th 



m e cw 2 K^T' 1 
2ne 2 nl ( K 1 (T I T 1 



with approximate multipliers 



cos 9 ■ g(X e ), 
+ 6T e 



sin 



h(X e 



(58) 



g(x e 

h(X e 



0.111n(l + 0.035X e 



2.011 exp 



X 



4.7 



exp 



xl- 2 
2.73 



0.011 exp (-||) 



(59) 



5 APPLICATIONS 

We computed the response tensor in uniformly magnetized relativistic plasmas with isotropic particle distributions. We found 
Faraday conversion, Faraday rotation, and absorption coefficients by numerical integration in the complex plane. Then we 
discussed the properties of natural modes of cyclo-synchrotron radiation and presented the results for specific electron energy 
distributions. We provided accurate practical fitting formulae for Faraday conversion and rotation. The method of complex 



plane integration allows to calculate both absorption and propagation coefficients consistently. In practice, formulae in lSazonov 



1 1969h are good enough for absorption coefficients. Therefore, we focus on improving the calculations of Faraday conversion 



and rotation coefficients. Faraday conversion and rotation coefficients can be generally expressed as pq oc h(~f, \B,6)n e B 2 \ 3 
and pv oc g(y, XB, 6)n e BX 2 , where h, g are functions of electron Lorentz factor 7, product of wavelength A and magnetic field 
B = \B\, and angle 8 be tween wavevector fc and magnetic field B. 

Previous work (e.g., Melrose! 1997bl ; Quataert fc Gruzinov 2O10h has shown that g w 1, if 7 « 1 for typical A, B, 6, but 



g <C 1 if 7 3> 1. This means the Faraday rotation becomes less important, when the electrons become relativistic. Therefore, 
the direction of LP plane changes little. Our calculations suggest the function g and the Faraday rotation coefficient were 
previously computed precisely, but function h and Faraday conversion coefficient were not. The function h grows from small 
to intermediate 7, but h steeply decreases, if 7 grows more. This means the Faraday conversion also becomes less important 
when the electrons become relativistic, although the peak Lorentz factor depends on frequency ratio and the particle 

distribution. In sum, as the absorption coefficients at a given frequency tj also decrease with 7, ultrarelativistic electrons 
interact less with radiation field. 

We have shown that Faraday conversion coefficient should be computed more precisely. Let us now demonstrate that 
imprecise estimates of propagation effects result in large error in polarized simulated fluxes for accretio n onto compact objects. 



We c ompare polarized spectra for two cases: accurate general relativistic polarized radiative transfer ( Shcherbakov fc Huanel 



2011) and assuming that Faraday conversion and rotation coefficients are given by linear approximations. We test two types 



of dynamical models for Sgr A* accretion to prove the case. Both model s assum e thermal electron distribution. 

First, we adopt the analytic model from iHuang. Takahashi fc Shen l|2009h . This model is established for Sgr A* based 



on radiatively inefficient accretion flow (RIAF) solution with different temperatures in ions and electrons. It reasonably fits 
Sgr A* polarized mm/sub-mm spectrum. On the left panel of Figure [9] we plot the accurate ratio Pq/pv (solid lines) and 
the ratio pQ,nn/ pv.Un in linear approximation (dashed lines) as functions of electron temperature T e for Q.o/uj = 10~ 4 and 



J Note, that the basis vectors in IShcherbako"^ <200Sh are different, thus a different sign of pg. 
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Figure 9. Left: Examples of accurate ratios of Pq/pv (solid) and their linear approximations pQ u n / Pv lin (dashed) for thermal electrons 
with Qo/lu = 10 -4 and different 9. Right: Simulated CP and LP fluxes based on accurate Fara day conversion/rotation (solid) and their 
linear approximations (dashed) for Sgr A* accretion model in (Huang, Takahashi & Shcn 2009) with inclination angle of 75° (black) and 
90° (red). 
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Figure 10. Simulated LP (on the left) and CP (on the right) fractions for Sgr A* for various prescriptions of Faraday rotation and 
Faraday conversion: for linear ap proximations pp i \ „ / py i \ „ (dashed red) and for accurate py and pQ (solid green). We employ the 
best-fitting dynamical model from lShcherbakov et all j2010T ) with dimensionless spin a* = 0.9. 



different k — B angles. On the right panel, we plot simulated circular and linear polarization spectra computed with accurate 
Faraday conversion and rotation coefficients (solid lines) and the results for the same dynamical model for assumed linear 
approximations of propagation coefficients (dashed lines) . The inclination angles of the disk i = 75° (black) and i = 90° (red) 
are chosen. 

The frequency ratio Oo/tu ~ 10 -4 corresponds to sub-millimeter band close to the event horizon of Sgr A*. There the 
electron temperature is T e ~ 10 10,5 ~ 11 ' 5 K. The Faraday conversion becomes strong (pq/pv > 1) at large angles 9 > 45°. 
One can also note from the right panel of Figure [4] that rjy absorption coefficient is less than 1% of rji. Therefore emissivity 
in V is equally weak according to Kirchhoff 's Law. So Faraday conversion plays a major role in the production of circular 
polarization. However, po/pv reaches a peak at a specific temperature around 10 n K then decreases again, while pQ,iin/ pv,Un 
monotonically increases to much greater than 1. Therefore, the amplitudes of circular polarization predicted for accurate 
(pQ,pv) are less than half of those predicted for approximate (pQ.im, pv.iin), although accurate and simplified propagation 
coefficients predict similar linear polarization. 



Changes in simulated polarized fluxes are shown on Figure[l0]for the best-fitting Sgr A* accretion model from lShcherbakov et al 



(2010). This model is inspired by three-dimensional general relativistic magneto hydrodynamic simulations. The structure of 
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magnetic field, velocity and density fields are taken directly from simulations. All free parameters are adjusted to achieve the 
best fit. Similarly to aforementioned analytic model, circular polarization in sub-millimeter band is significantly lower, when 
the precise Faraday conversion and rotation coefficients are adopted for the same dynamical model. Thus substantially lower 
predicted CP fractions, when the precise propagation coefficients are adopted, is a generic model-independent result. Our 
calculation showed that linear approximations of Faraday conversion/rotation are invalid for electrons with such high energy. 
They significantly overestimate the circular polarization for the relevant range of ratios fio/w. When fitting polarized obser- 
vations and predicting polarized spectra one should adopt accurate Faraday conversion and rotation coefficients computed in 
§ 12.21 or the simplified fitting formulae provided in § 14.21 

The validity of linear approximations of Faraday conversion and rotation coefficients depends not only on the energy 
of electrons, but also on the frequency ratio Qo/ui. The observational frequency in near IR (NIR) is ~ 10 3 larger than that 
in sub-mm. The typical frequency ratio fio/w is around 10~ 6 — 10~~ 8 , which yields the values of (pQ,pv) similar to those 
°f (PQ.iin, Pv,iin) for 7 < 50. Thus, the linear approximations of Faraday conversion and rotation coefficients still adequately 
describe the correspondent effects in NIR. 

Observations of sub-mm circular polarization from Sgr A* are one of tod ay's big interests a nd challenges. CP fraction 
is only about 1% in sub-mm, but was alr eady detected at several frequencies IjMunoz et al.ll201lr ). Still, there is substantial 
spread between different Sgr A* models 1 Shcherbakov et al. 2010) in CP fluxes at frequencies, where CP flux was not yet 
measured, for example 88 GHz, f45 GHz, and 690 GHz. Observations at these frequencies can help to further discriminate 
between models of various types, which have different black hole spins. 
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Appendix A: Definitions and derivations 



Here we summarize all definitions and derivations related to the response tensor. They are similar to those in ( Melrose! 

1997al lbf), except for the signature of the metric tensor and the geometry of the coordinate system. We list all related tensors 
which are different from those in Melrose's work. 

Metric tensor in Minkowski space-time: g^" = diag[— 1, f, f, 1], 
Wave vector: fc M = w(l, 0, 0, 1) T . 
Tensor of magnetostatic field: 



B 



1/2 



( 


— cos ( 

\ -suit 



Auxiliary tensor for electron velocity perturbations: 
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Auxiliary tensor for electron position perturbations: 
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In 6-dimensional phase space: electron velocity v, electron momentum m e p, and distribution function of electrons f(p). 
In 8-dimensional phase space: Lorentz factor 7 = V\/ 1 — |v| 2 , velocity of electron W = (7,u) = (7,7V), momentum 
7-" 1 = (jm e , m e p) = ( r ym e , jm e v), and distribution function of electrons T(P) = 2m e 8 (V 2 + m 2 )f(p) . 
The general form of the response tensor for an arbitrary isotropic distribution of electrons is 
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Appendix B: Derivation of the response tensor for monoenergetic particle distribution 

Let us introduce W = [1, 0, 0, 0] — 4-velocity in the rest frame of plasma. Then 
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which coincides with Eq. [9] We denoted 7?. M (£) = k\T^ x {S > ) and introduced an auxiliary variable 5 M . 
Appendix C: Derivation of the response tensor for thermal particle distribution 

Substituting the thermal distribution into Eq. [12] we obtain 
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where 

IpI = \Ji 2 - i- We can simplify the above formula as follows: 
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Following Eq. [9] and using 
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which is Eq. l57l 
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